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Diagrammatic theory for Periodic Anderson Model has been developed, supposing the Coulomb 
repulsion of /— localized electrons as a main parameter of the theory. /— electrons are strongly 
correlated and c— conduction electrons are uncorrelated. Correlation function for /— and mass 
operator for c— electrons are determined. The Dyson equation for c— and Dyson-type equation for 
/— electrons are formulated for their propagators. The skeleton diagrams are denned for correlation 
function and thermodynamic functional. The stationary property of renormalized thermodynamic 
potential about the variation of the mass operator is established. The result is appropriate as for 
normal and as for superconducting state of the system. 

PACS numbers: 71.27.+a, 71.10.Fd 



I. INTRODUCTION 

The study of the systems with strongly correlated 
electrons has become in the last time one of the cen- 
tral problem of condensed matter physics. One of the 
most important models of strongly correlated electrons 
is periodic Anderson model (PAM)W. This model is 
used to describe the physics of mixed valence systems, 
heavy fermion compounds, high-temperature supercon- 
ductivity as well as other phenomena in which the strong 
Coulomb repulsion of the localized electrons is present. 
This model describes the intermetallic compounds which 
contain magnetic moments of rare earth or actinide ions 
included in the host metal. This ions have a partially 
filled /— shell and can be considered as scattering cen- 
ters for conduction electrons of the host metal. Because 
of the strong Coulomb repulsion of the electrons with op- 
posite spins located at the same site of lattice the mag- 
netic ion electrons are strongly correlated. There is also 
the hybridization of states between the uncorrelated con- 
duction electrons and localized correlated ones when both 
of them are present on the same lattice site. Magnetic 
properties of the impurities ions affect in a different man- 
ner the properties of the host matrix and of the system 
as a whole. For different regime of physical parameters, 
determined by Coulomb local interaction, hybridization 
of the wave functions and exchange interaction, it is pos- 
sible to obtain different classes of the system phases. 

There are already an enormous number of approximate 
methods and approaches devoted to PAM, as perturba- 
tion expansions, static and dynamic mean field theories, 
variational and numerical approaches, large N expan- 
sion , slave boson methods, non crossing approximations 
(NCA), Bogoliubov inequality method and others. Also 
some exact results are known, obtained in special with 
the Bethe ansatz, renormalization group methods and 
Bogoliubov inequality method. We will not enlarge upon 



the most essential stages in the development of this model 
because exists a number of consistent reviews I 2-7 ' and 
books I 8,9 ' on this field and we shall use the references to 
previous our papers. 



II. MODEL HAMILTONIAN 

We consider the simplest form of PAM with a spin de- 
generation of the level of localized /— electrons, a simple 
energy band of conducting c— electrons, Coulomb one- 
site repulsion U of correlated /— electrons with opposite 
spins and one-site hybridization between both group of 
electrons of this system. The hamiltonian of the system 
reads: 
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Here V is the hybridization amplitude assumed con- 
stant. We have indicated with Cf a (fJ~ a ) the creation op- 
erator for an uncorrelated (correlated) electron with spin 
a and i lattice site, m a is the number operator for /— 
electrons, e(k) is the band energy with momentum k of 
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conductivity electrons spread on the entire width W of 
the band, e/ is the energy of localized electrons. Both 
these energies are evaluated with respect to the chemical 
potential \i. 

The approach proposed in this paper generalizes the 
diagrammatic theory of normal and superconducting 
phases of strongly correlated systems proposed in pre- 
vious papers I 10-19 !. 

The strong on-site repulsion U between /— electrons 
of opposite spins is the main term in the Hamiltonian. 

As the conduction electrons can belong not only to the 
s— but also to the d— atomic shell, their Coulomb re- 
pulsion can also be important. In this case the extended 
PAM must be used [ 2 °1. For simplicity the correlations 
of c— electrons are not considered and one subsystem 
is of c— uncorrelated and the second of /— correlated 
electrons. Because of strong localization of the /— elec- 
trons they cannot hope from one lattice site to another 
and their derealization is due only to the hybridization 
of the /— and c— states with matrix element V. It is 
obvious that at V 7^ in the given model with two sub- 
systems superconductivity arises simultaneously in both 
subsystems. 

In the present paper we develop the thermodynamic 
perturbation theory for the system in the superconduct- 
ing state with Hamiltonian (1) under the assumption that 
the term responsible for hybridization of c— and /— elec- 
trons is a perturbation. 

The Hamiltonian H® of the uncorrelated c— electrons 
is diagonal in band representation, where as the Hamil- 
tonian Hj is diagonalized by using Hubbard transfer op- 
erators I 21 ' . Therefore in the zeroth-order of the pertur- 
bation theory the statistical operator of grand canonical 
ensemble of the system is factorized in the momentum 
representation for c— and in local representation for /— 
electrons: 

exp[-(3(H° c +H ( >)} = [JexpH^k)) x 
x JJexp(-^(i)), 

i 

H° c Jk) = e(k) C+ C kff , (3) 

cr 

We use the series expansion for evolution operator: 

P 

U{(3) = TcM-J H mt (r)dT). (4) 


in the interaction representation for electron operators 
(a = c, /): 

a(x) = e TH " a(x)e- TH ° ,a(x) = e rH V(x)e- Tff °. (5) 

Here by x means (x, a, r). 

We shall denote by (TAB...) Q the thermodynamic av- 
erage with zeroth-order statistical operator (3) of the 



chronological product of electron operators (AB...). Such 
averages are calculated independently for c— and /— op- 
erators with using for c— electrons the Wick Theorem of 
weak quantum field theory and by using for /— electrons 
the Generalized Wick Theorem (GWT) proposed by us 
in papers [ 10 ~ 13 1 for strongly correlated electron systems. 

In the superconducting state, unlike the normal one, 
nontrivial statistical averages of operator products with 
even total number but inequal number of creation and 
annihilation electron operators are possible. They re- 
alize the Bogoliubov quasi-averages I 22 ' or Gor'kov I 23 l 
anomalous Green's functions. To unify the calculation 
of statistical averages for normal and superconducting 
phases it is useful to assign an additional quantum num- 
ber a, called by us charge number I 15 !, with the values 
±1, which can be add to electron operators according the 
rule (a — c, /): 

a«(x) = { a = l \ (6) 

v ; [ a + (x), a = — 1. v ' 

In this representation the interaction operator H int be- 
comes: 

Hint = v5>/r°C£. (7) 

Obviously, introducing a new quantum charge number 
leads to additional summation over it values in all dia- 
gram lines and to an additional factor a in the vertices 
of diagrams. 

Now, after such introducing, it is irrelevant whether 
one deals with creation or annihilation operators. First 
of all we shall enumerate the main results of diagram- 
matic theory obtained in the previous paper I 15 ! neces- 
sary to our proving of stationary theorem. Such theorem 
for uncorrelated many-electron systems in normal state 
has been proved by Luttinger and Word I 24 ' . 

III. PERTURB ATIVE TREATMENT [15] 

We use the definition of the one-particle Matsubara 
Green's functions for c— and /— electrons 

<w(s|*') = - (Tc a (x)c- a >(x')U(f3)) c , 

(8) 

gH*W) = (T.r(x)r a '(x')u(p));, 

where index c for means the connected of the di- 

agrams which are taken into account in the right-hand 
part of definition (8). 

The following condition is fulfilled 

G*A*W) = -G a _ a>! _ a (x'\x),a = (c,f). (9) 

Between this new definition and traditional one I 23 ' 
there is a relation 

G a hl {x\x') = G a (x\x'), 
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Gi-i{x\x') = F a (x\x'), 
G a _ M (x|x') = F a (x\x>), 
Gl^ixlx') = -G a (x'\x). (10) 

In the presence of strong correlations of /— electrons the 
(GWT) contains additional terms namely the irreducible 
one-site many-particle Green's functions or Kubo cumu- 
lants of the form (x = x, er, r) : 

G^ ir [a 1 ,x 1 ;...;a 2n ,x 2n } = (Tf^.-.f^ = (11) 
= <W..£x lX2n (T/^(r 1 ).../^;(r 2n ); r , 

where in simplest two-particle case we have the following 
definition of the irreducible function 

<T/^(r 1 )/ ff Q 2 2 (r 2 )/^(r 3 )/ ff Q ;(r 4 )); r = (12) 

= (T/-(r 1 )/-(r 2 )/-(r 3 )/-(r 4 )) - 

-[(T.Cl(T 1 )f^(T 2 )) (r.^(r 3 )C(r 4 )) + 

+ (TC (n)C (r 4 )) (r 2 )/£ (r 3 )) - 

- (T/:(ri)C(r 3 )) (^(T 2 )C(r 4 ))J. 

For n > 3 the irreducible Green's function contains in 
its right-hand part besides the products of one-particle 
propagators also their products with irreducible functions 
of smaller number of particles. There are present also the 
product of irreducible functions Gn} lr and G^° 2 lr with 
condition n\ + n 2 = n and so on I 10 ^ 13 !. 

As a result of applying these theorem we obtain for the 
rcnormalizcd conduction electron propagator the contri- 
butions depicted on the Fig. 1 The contributions of the 
diagrams Fig. 1 b) and e) are the following 

V 2 E E^^Gl°)(-|l)G^ 2 (l|2)G^( 2 l-'), 

a\d2 12 

^ £ ^a ia2 a 3 «4Gl )(x|l)G^) 3 (2|3)x 

ai...a4 1...4 

x<°>N') (Tf^f 2 a2 f^U ai )l r 

correspondingly. It demonstrates the dependence of the 
diagrams from the charge quantum number a. 

The contributions of perturbation theory for /— elec- 
tron propagator are depicted on the Fig. 2 The contri- 
butions of diagrams Fig. 2 b) and c) are equal to 

V 2 E E«i^G^ 1 ) (x|l)G^) 2 (l|2)GS(2|x'), 

C1102 12 

— J2 E a i a 2G^) 2 (i|2)(r/«/r Qi /r.c Q ') o , 

QlQ2 12 

correspondingly. 

Between the diagrams for the one-particle /— propaga- 
tor there are strong and weak connected ones. The weak 
connected diagrams can be separated in two parts by cut- 
ting one propagator line. The sum of all strong connected 



diagrams for /— electron belong to the correlation func- 
tion which is denoted by us as A aa >(x\x') function. The 
quantity A aa /(a;|a;') is defined by the equation 

A aa >(x\x>) = G f J°)(x\x') + Z aa ,(x\x'), (13) 

where the function Z aa i(x\x') contains the contribution 
of strongly connected diagram based on the irreducible 
many-particle Green's functions. 

The strong connected part of the c— electron propa- 
gator without the external lines is determined by us as a 
mass operator for uncorrelated electrons. This quantity 
is denoted as Y, aa i(x\x'). 

A simple relation exists between these two functions: 

£ QQ ' (x\x r ) = V 2 aa'A aa , (x\x r ) (14) 

The analysis of the propagator diagrams permits us to 
formulate the following Dyson equation for uncorrelated 
electron propagator 

G c aa ,(x\x') = G c a {0 J(x\x')+ (15) 
+ E E G ^^l 1 )^ 1 « 2 ( 1 l 2 )^ 2Q '(2N'). 

aia2 12 

At the same time we can formulate the Dyson-type equa- 
tion for correlated electron propagator G^ : 

G f aa ,(x\x') = A aa ,(x\x') + (16) 
+ V2 E E a i^A QQl (x|l)G^2 2 (l|2)G^ Q ,(2| a; '). 

<X\Ct2 12 

In equation (15) and (16) as in the previous equations 
which contain repeated indices (1,2) is supposed sum- 
mation by sites indices (1,2), spin indices (<7i,<7 2 ) and 
integration by time variables (n , t 2 ) in the interval (0, /3). 

Unfortunate the Dyson-type equations far correlation 
function A aa > and mass operator £ QCe < don't exist. There- 
fore the calculation of the c— and /— renormalizcd propa- 
gators needs the approximations based on the summation 
of special classes of diagrams. 

On the Fig. 3 the skeleton diagrams for the correla- 
tion function A aa >(x\x') are depicted . They demonstrate 
impossibility to formulate Dyson-type equation for this 
function. The number of skeleton diagrams depicted on 
the Fig. 3 for the correlation function is infinite. 

The contribution of the diagrams Fig. 3 b) and c) is 
the following: 

\ E E«i a ^ lQ2 (i|2) (Tf2fr ai f? 2 f-, a y r o , 

a±a2 12 

~\ E E a i a 2a 3 a4G^ lQ2 (l|2)G^ Q4 (3|4)x 

ai...«4 1...4 

x (T/ x Q /r ai / 2 a2 / 3 " a3 /4 a4 /,- Q ')^ , 

correspondingly. 
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FIG. 1: The first six orders of perturbation theory for conduction electron propagator. The solid and dashed thin lines depict 
zero order propagators for c— and /— electrons correspondingly. The rectangles depict the irreducible Green's functions. The 
points of diagram are the vertices with a and V contributions. 
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FIG. 2: The contributions of the first four orders of perturbation theory for the /— electron propagator. 



If we take into account only the first term of the right- 
hand part of Fig. 3 we obtain the simplest Hubbard I 
approximation with consideration only of the chain-type 
diagrams. 

The diagram Fig. 3 b) is the simplest contribution 
to correlation function which takes into account the elec- 
tronic correlations. The diagrams Fig. 3 b), c) and d) are 
localized and their Fourier representations in real space 
are independent of momentum. There are also other dia- 
grams of this kind with irreducible functions G 5 , G 6 
and so on. The coefficients before these diagrams are de- 
termined by the number t^t^t , where 2n is the perturba- 
tion theory order of diagram. The last diagram of Fig. 
3 is not local and its Fourier representation depends of 
the momentum. To dynamical mean field theory only the 
first group of local diagrams of Fig. 3 correspond. 

The transition of the diagram contribution from su- 
perconducting version to the normal one is realized by 



the condition of equality to zero of the sums of all a— 
indices of every dynamical quantity. For example such 
transition of the diagram Fig. 3 b) is conditioned by the 
equalities a± — ct2 — and a — a\ + cti — a 1 ~ with 
solution a\ — 0:2 and a = a?. 

The summation by u\ gives us two equal contributions 
and the coefficient before diagram increases twofold and 
becomes I instead originally ^. 

In normal state the correlation function A(x|a;') has a 
form depicted on the Fig. 4. 

The new coefficient before the diagrams take into ac- 
count the existence of different possibilities of transition 
from superconducting to normal state. 

After discussion of the propagators properties we shall 
proceed to the main part of our paper and investigate 
the properties of evolution operator average. 
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FIG. 3: The skeleton diagrams for correlation function A aa i(x\x'). The thin dashed line is zero-order /— electron Green's 
function. The rectangles depict the many-particles irreducible Green's function. The double solid lines depict the renormalized 
conduction electron Green's function G c aa ,(x\x'). 



A(x | x') = 




FIG. 4: Correlation function A(x\x') in normal state. All the lines correspond to normal propagators and have a direction of 
propagators and of the arrays in the vertices. 



IV. VACUUM DIAGRAMS 

By using the perturbation theory we have obtained 
for the connected part of evolution operator average the 
contributions depicted in the Fig. 5. 

The contributions of the first, third and fourth dia- 
grams of Fig. 5 are enumerated below: 

" T ^ E^ a ^ 1 l(l|2)G^°i 1 (2|l), 

Q1Q2 12 

X £ E a i a 2a 3 a 4 G^) 2 (l|2)G^) 4 (3|4)x 

Qi...a4 1...4 

x(r/r7 2 "7 3 "7" 4 ); r , 

E E a i a 2a3a4a 5 a 6 G^ 2 (l|2)G^) 4 (3|4) x 

ai...a 6 1...6 

xG^im (Tfr>fpf<i a3 ftfr 5 f6 6 )o • 



' 48 



In normal state the diagrams of Fig. 5 are changed. 
The direction of propagator lines and of the arrays at 
vertices points appear together with new coefficients be- 
fore the diagrams. This changing is demonstrated on the 
Fig. 6, where some of the diagrams of Fig. 5 are demon- 
strated. Vacuum diagrams in superconducting and nor- 
mal states contain the factor — , where n is the order 
of perturbation theory in which given diagram appears. 
This factor makes difficult the investigation of this contri- 
butions. In order to remove this coefficient it is necessary 
to use the trick of integration by constant of interaction 
V. The result of such integration is depicted on the 
Fig. 7. 

On the base of series expansions for renormalized prop- 
agators of the conduction c— electrons (see Fig. 1), lo- 
calized /— electrons (see Fig. 2) and definition of the 
correlation function A aa ' (x | %') we can prove that the 
contribution of the integrant of Fig. 7 in every order of 
perturbation theory can be presented itself as the prod- 
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FIG. 5: Vacuum diagrams of first eight orders of perturbation theory in superconducting state. 
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FIG. 6: Some of the vacuum diagrams in normal state of the system. 



uct of some contribution from G c and some one from A. 
If the contribution of G c is of m order of perturbation 
theory and of A of n-i order when the order of (C/(/3))n 
is equal to n with the condition n\ + n% + 2 = n which 
must be satisfied. There are different possibilities to sat- 



isfy this condition and all of them must be taken into 
account. 

For example there are three possibilities to compose 
from G c and A the sixth diagrams of Fig. 7. These pos- 
sibilities are enumerated below on the Fig. 8. Other ex- 
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FIG. 7: Vacuum diagrams of the system in superconducting state after integrating by interaction constant. 



ample of vacuum diagram of eighth order of perturbation 
theory is presented on the Fig. 9. Only all the three pos- 
sibilities give us the correct coefficient — -g- — i — -i = — i of 
diagram. Other possibilities don't exist. These examples 



demonstrate the general statement that the integrand of 
the evolution operator average can be presented itself as 
a product of A 2 G C A of the form 



V 



Y J Y J aa ' x2GC ^'( x \ x, \ x ) Aa ' a ( x '\ x \ x ) 



(17) 



V V 

j dXj2J2G C aa^W\m»'a(x'\x\X) = - J ^Tr(G c (A)£ c (A)). 



where the operators G C (A) and S C (A) have the matrix 
elements G c aa , (x\x' \X) and Sac/ (a; | a:' | A) correspondingly. 
Index A underline that these quantities depend of the 
auxiliary constant of integration A. 

Therefore the thermodynamic potential of our system 
F is equal to 

v 

F = F + - p J ^Tr[G c (A)£ c (A)]. (18) 
o 



This expression for renormalized thermodynamic poten- 
tial of the strongly correlated system contains additional 
integration over the integration strength A and is awk- 
ward because it. Equation (18) generalizes the result of 
Luttinger and Ward I 24 ' proved for non-correlated many- 
electron system in normal state. 

Our generalization has been obtained for the case of 
strong correlations of special kind which contains one un- 
corrected subsystem and one strongly correlated and we 
admit also the existence of superconductivity in both of 
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FIG. 8: Three possibilities to organize the vacuum diagram of sixth order of perturbation theory. The correct coefficient § is 
obtained by summing all the possibilities. 



them. 

Luttinger and Ward have proved the possibility to 
transform this expression into much more convenient for- 
mula without such integration. For that they used a 
special functional constricted from skeleton diagrams the 
lines of which are the renormalized electron Green's func- 



tions. We shall use the skeleton diagrams of strongly 
correlated system which differ essential from Luttinger 
and Ward I 24 l case and transform equation (18) to more 
convenient form. 

In our strong correlated case we introduce the following 
functional 



Y = -^-Tr{\n(G c{a) t - 1) + G c £} + Y', 



(19) 



r 



which is the generalization of the Luttinger- Ward ^ 24 1 
equation just for the strongly correlated systems. Here 
operation Tr use the summation by a, a, i and integra- 
tion by r. 

The quantity Y' contains all peculiarities of the 
strongly correlated systems and is presented itself as a 
sum of closed linked skeleton diagrams, constructed from 
irreducible Green's functions of correlated electrons and 
full Green's functions of uncorrelated electrons. 

On the Fig. 10 are depicted some of simplest skeleton 
diagrams for functional Y'. These diagrams depend of 
the interaction strength V not only through the factors 
in the front of each diagram but also through the depen- 



dence of full Green's function G C (V). The contributions 
of the diagrams a), b) and c) are following 
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FIG. 9: Three possibilities to obtain one of the vacuum diagram of eight order of perturbation theory. 
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FIG. 10: Skeleton diagrams for functional Y' . 
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From Fig. 3 and Fig. 10 it is possible to demonstrate 
the following equation: 

8j3Y' Q\Ol2 2 a / s 1 

( 2 °) 

If we take into account only the explicit dependence of the 
functional Y' of interaction constant V without consid- 
ering the dependence of the full Green's functions G C (V) 



,dY' 



aia2 12 



(21) 

Because of the necessity to have the functional derivatives 
over mass operator E we shall use the Dyson equation 
(15) rewritten in the form [a; = (a, x, a, t)]: 

G c(0)_1 (x,x') = G°~ x {x,x') + T,(x,x') (22) 
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We obtain 



and 



6G c (x,x') 

snv,y') 



G c (x,y)G c (y',x'), 



(23) 



-Tr\n(l-G c ^±) = -G c (y',y), 



-Tr(G c t) = G c {y',y) + {G c ±G c ) y , y (24) 



S^(y,y 

By summing these equations we obtain: 

^(4 )Tr{ln( ^ (0)£ - 1) + ^ } = 

= -±{&±&) y , v . (25) 

On the base of definition of the functional Y' Fig. 10 and 
equation (20) we have 

SY' _ y SY' SG c (x,x') _ 1 , c -, 
tfEfo.y') ^6G'(x,x>) 5^{y lV ') 2(3 [ h ' v ' 

(26) 

As a result we obtain the stationarity property of the 
functional Y: 



6Y 



6V(y,y>) 



0. 



(27) 



Now we shall discuss the derivative over interaction 
constant V of functional Y. We shall taken into account 
the stationarity Y about E and G c and equation (21). 
We obtain: 

y— = | E = v 9 -^ | E = ^1(28) 

|E .sew ay |E /3 1 ; 



From equation (18) we have: 

dF_ _ Tr{G c ±) 
dV ~ ' 

and as a consequence we establish 
dV dV 

with the solution 

F = Y + CONST. 

This constant is Fq. Therefore 

F = F + Y, 
with the stationary property 
5F 



0, 



SH(x, x') 

as in superconducting and in normal states. 



(29) 

(30) 

(31) 
(32) 

(33) 



V. HEAT CAPACITY [25] 

As a illustration of our results we shall consider the 
problem of finding the heat capacity of our strongly cor- 
related system in normal state and at the low tempera- 
tures. 

The heat capacity at constant volume V is equal to 



c 'v - T(—) v , 
where the entropy S is given by 



(34) 



(35) 



The quantity \i is the chemical potential at temperature 
T. 

At low temperature we may expand F(fj,, V,T) and // 
in even powers of the temperature: 

F = F(^ o ,V,0) - ^(fi ,V)T 2 + 

M(T) - Mo+M'T 2 + .., 

where /xo is the value of chemical potential for correlated 
system at T = 0. 
As a result we have 



C v = l(ji ,V)T + ... 



(36) 



the linear dependence of the heat capacity of the tem- 
perature at low its values. Therefore to evaluate the co- 
efficient j(fj, ,V) it is necessary to obtain the expansion 
of F in powers of T by using the expression (31), (19) 
and Fig. 10 for Y' obtained by us in previous part of the 
paper. 

We know that the expression (19) for Y and (31) for 
thermodynamic potential is stationary with respect to 
changes in the proper self-energy Ek(iwo)- Therefore be- 
cause we are interested in the first corrections to F, we 
can neglect the explicit temperature dependence of mass 
operator E(k|iw) and G c (k\iuj) and replace them by val- 
ues E(k|iwo)|T=o and G c (k|iwo)|T=o calculated at T = 0. 
Thus the first correction to the T = value of F (31) 
comes only from the difference between the u n sums in 
expression (31) and what we would get if we replace them 
by integrals according to the equation 



-Y = - 

2tt 



did. 



(37) 



Now it is necessary to consider the functional Y' on 
Fig. 10. Since each line of a skeleton diagram of func- 
tional Y' contains an uj n sum, the total first correction 
to Y' is obtained by correcting the computation in each 
diagram for a single line and use equation (37) for the 
other uj n sums of the diagram, finally summing over ev- 
ery line. 
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In such a way we obtain the contribution of one line 
of skeleton diagram multiplied by the number of lines of 
skeleton diagrams. This number changes the coefficient 
before the skeleton diagram of Y 1 and new coefficients 
correspond to new contribution to Y' equal to the self- 



We use the Poisson equation for the uj n sums and write 



where C is contour which surrounds in anti clock wise di- 
rection the poles of the function (eP z + 1) _1 in the points 
z = iuj n = ( 2 "+ 1 ) 7r i. The term in last equation propor- 
tional to T 2 is obtained by usual Sommerfeld technique. 

The details of such computation will be discussed in 
other place. 



VI. CONCLUSIONS 

We have developed the diagrammatic theory for PAM 
on the base of new conceptions proposed by us for 
strongly correlated electron systems. 

We introduced the notion of correlation function 
A aa r(x\x') of /— electrons (see Fig. 3) which is the in- 
finite sum of strong connected irreducible Green's func- 
tions and which contains the most important spin, charge 
and pairing fluctuations of the correlated /— electrons. 
This correlation function determines the mass operator 
{%\%') (14) of the uncorrelated conduction electrons. 
The both these quantities A and X permit us to formu- 
late the Dyson equation for c— electrons (15) and Dyson- 
type equation (16) for /— electrons. These results are 
expressed in general form appropriate as for normal and 
as for superconducting state. 

We have obtained the skeleton diagrams for A func- 
tion and demonstrated their dependence from irreducible 
many-particle Green's functions Gl°' 4r (l, In) with all 
values of n and also of c— electrons full propagators. 
Thanks the presence of these irreducible Green's func- 



energy one. This quantity has the form of second term in 
right hand part of functional (19) for Y having the oppo- 
site sign. When we combine the both part of functional 
Y these quantities are reciprocally canceled. 

Finally we obtain to the first order for F the equation 



(38) 



them as an integral 



(39) 

I 

tions it is impossible to formulate Dyson-type equations 
for A and X quantities. 

The results are appropriate as for normal and as the su- 
perconducting state. Unification of the investigation for 
the both phases was possible thanks the introducing of 
the notion of quantum charge number a and the rewrit- 
ten of the interaction Hamiltonian in such new form. 

From Fig. 3 it is clear that the simplest contribution 
that takes into account /— electron correlations is re- 
duced to first two terms of right-hand part of this figure. 
All the terms of Fig. 3 besides the last one and also 
other omitted diagrams like them are local with Fourier 
representation independent of momentum. These terms 
correspond to the structure of dynamical mean field the- 
ory. Last diagram of Fig. 3 and other more complicated 
diagrams with more number of irreducible Green's func- 
tions depend of momentum and take in consideration of 
the space fluctuations. The local contributions take into 
account only of the fluctuations in time. 

We have demonstrated the transition of our diagram 
from superconducting to normal state by using the addi- 
tional conditions imposed on the charge quantum num- 
bers of which depend the dynamical quantities. 

The special investigation of vacuum diagram has been 
done after introducing the auxiliary interaction strength 
and integration by it of these diagram contributions. We 
have proved that this integrant is equal to the product 
of two matrices T, aa > and G c aa , . 

Then we have introduced special functional in the form 
of skeleton diagrams and proved it coincidence with ther- 
modynamical potential. This expression has the property 



J 



F = F - — ^^ln[G c (°)^„)S(k| iW „)| T= o - l]e^" 0+ = 

£rk U> n 

= ^o + ^EE ln ^« - - ^ E E l*[£(k|ia; n )|r=o " ( GC( ° } (*"»))" V 

trk u» trk Ul n 

I 



J 



F = F o - 2 E _/ ^T^[G c(0) WS(k|z)| T=0 -l], 
(Tk c 
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of stationary relative the changing of the mass operator 
or full Green's function G c of conduction electrons. 
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